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Abstract 

We consider one-loop scalar and tensor integrals with an arbitrary number of external legs 
relevant for multi-parton processes in massless theories. We present a procedure to reduce Ap- 
point scalar functions with generic 4-dimensional external momenta to box integrals in (4 — 2e) 
dimensions. We derive a formula valid for arbitrary N and give an explicit expression for N = 6. 

Further a tensor reduction method for iV-point tensor integrals is presented. We prove 
that generically higher dimensional integrals contribute only to order e for N > 5. The tensor 
reduction can be solved iteratively such that any tensor integral is expressible in terms of scalar 
integrals. Explicit formulas are given up to N — 6. 
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1 Introduction 



At future colliders multi-particle /jet final states will become more and more dominant. A precise 
theoretical description of the QCD reactions is desirable in order to have a better control on the 
backgrounds for various search experiments. Especially at hadron colliders lots of multi-jet data will 
be collected which have to be confronted with theoretical predictions. This motivates the accurate 
calculation of multi-parton reactions. 

For tree level calculations the construction of matrix elements with large numbers of particles 
in the final state is well established However, tree level results are very unstable with respect 
to variations of the renormalization and factorization scales. As 2 — > N parton cross sections 
behave as otg, at least next-to-leading order precision is necessary to stabilize the predictions. The 
knowledge of the respective one-loop corrections is thus mandatory. For 2 — > 3 parton scattering 
NLO contributions have been calculated in |@, [|f|], for e + e~ — * A jets see || ||. 

The main technical difficulties for constructing amplitudes consist in the treatment of the occur- 
ring TV-point scalar and tensor integrals. Because of infrared (IR) divergencies the standard methods 
°f @> H> @ are no ^ directly applicable. The authors of H used a formalism working in Feynman 
parameter space. The generation of parameter integrals with nontrivial numerators connected to 
tensor integrals was done by differentiation techniques. The formalism produced formally higher 
dimensional integrals which in the end canceled out. This cancellation had to be shown by explicit 
calculation. In this paper we present a proof that this cancellation mechanism is true for arbitrary 
N. Furthermore, a scalar reduction formula was derived in H which relates scalar integrals in dif- 
ferent space-time dimensions to each other. For 4-dimensional external momenta the formula could 
not be shown to hold true in general. The problem stemmed from the presence of singular matri- 
ces. We rederive the formula and treat the necessary inversions of these matrices by using so-called 
pseudo-inverse matrices. To keep the external momenta in 4 dimensions is advantageous because it 
allows to use helicity techniques . 

Another formalism based on the generation of tensor integrals by differentiation of scalar integrals 
serving as generating functionals has been proposed in jfl], [l2| . The reduction of the scalar integrals 
works basically in the same way as in |§ with the same problems in the case of 4-dimensional 
kinematics. Their method has been designed for massive integrals and the applicability to the 
massless case has not been worked out. The authors state that a regulator mass is needed for 
infrared divergencies in the case of 4-dimcnsional kinematics. 

A generalization of the Passarino-Veltman techniques dealing with the problem of vanishing 
Gram determinants has been developed in |T^| . The authors show explicitly how to reduce box 
tensor integrals to scalar integrals. Another approach which uses helicity methods for the reduction 
of tensor integrals to scalar integrals has been discussed in |l4], |l5| . 

The complete reduction to scalar integrals is also possible in our tensor reduction scheme. Our 
formalism is in a sense the generalization of the Passarino-Veltman methods used for the calculation 
of electroweak radiative corrections to the massless case. 

The paper is organized as follows. In section 2 we rederive a reduction formula for scalar integrals. 
We concentrate on the case of 4-dimensional kinematics from the start and show that any Ap- 
point scalar integral with N > 6 is a linear combination of pentagon integrals which in turn are 
combinations of box integrals plus terms of 0(e). We give an explicit formula for the 6-point function 
with all external legs on-shell. In section 3 we derive our tensor reduction formalism which combines 
Passarino-Veltman-like methods with Feynman parameter space techniques. We prove that for 
generic 4-dimensional kinematics all higher dimensional Appoint functions drop out for arbitrary 
N > 5 and generalize the reduction methods of || [l7| to arbitrary N. We construct a hierarchy 
of tensor formulas up to N = 6 and rank < N which can be solved by iteration. The explicit 
expressions and the derivation of some basic formulas are given in the appendix. 
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2 Reduction formula for massless scalar integrals 



In this section we will first derive a scalar reduction formula valid for an arbitrary number of external 
legs. Then we will use these formulas to derive explicit representations for scalar integrals up to 
N = 6. 

2.1 Derivation 

Consider a one-loop, scalar iV-point function with massless propagators for N > 4. If all legs are 
off-shell the integral is finite and can be treated in four dimensions B. If at least one external 
leg is massless it is infrared divergent and needs a regulator. In the framework of dimensional 
regularization, 4-dimensional methods are not applicable anymore. We work in n — A— 2e dimensions 
in the following with the external momenta kept in four dimensions. 

With the momentum flows as indicated in Fig. 0, we define the propagator momenta as qi = k—ri 
with ri = pi +r;_i for I from 1 to N and ro = tn- Momentum conservation allows to choose tn = 0. 




Figure 1: N -point graph. 
The corresponding analytic expression in momentum and Feynman parameter space is 



I N (R) = / dn 



1 



= (-1) N T(N -n/2) 



TTN 2 

IL=i if 



*(i-E£i*i) 



d N z 



(1) 



Herein R = (n,...,r N ), dn = d n k/(in n / 2 ). 
S which is related to the Gram matrix G by 

1 

~2 

2 r k ■ ri 



{z-S-z) N - n / 2 
The kinematic information is contained in the matrix 



Ski 
Gki 



1 



(n-f/t) = -^i G ki - vi - v k ) 



v k =G kk /2, M = l, 



(2) 



Although it is well known |l2|, [T^] that the A^-point integral can be split into a finite, (6 — 2e)- 
dimcnsional integral and a part with less external legs containing the infrared poles, we want to 
present our derivation which later allows to deal with the problem of vanishing Gram determinants 
in a transparent manner. As an ansatz we write (0) as a sum of (one-propagator) reduced diagrams 
and a remainder. 
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We want to show in the following that one can find coefficients bi such that Ifi n contains no IR 



poles. In standard Feynman parametrization one gets with D = E^li z i it 



ifin 



Y{N) j o d N z5(l-J2zi) J 



dn 



D 



N 



and after a shift k = & + E^Ii z z r ' 



JV 



I fin = T(N) / d N z6(l-J2zi) dk 



1=1 



(k 2 — M 2 ) N 
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z, q 3 



ij - zi)n 
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Now the term in square brackets in Eq. (ph can be written as 



1$ 



JV 



N 



= ~(fc 2 + M 2 ) bj + z j (l + 2(5 ' 

3=1 2=1 



If now the equations 

(5 • % 

are fulfilled, the second term on the right-hand-side of (|7|) vanishes and one finds 

k 2 + M 2 



-\, j = l,..;N 



N 



N 



1 = 1 



Ifin = -T(N) \J2 b n J Q d N z6(l-J2z l )jdkj: 
Finally the loop momentum integration gives 

Ifin 
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(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



The (6 — 2e)-dimensional integral is IR finite, as can be seen by a power-counting argument in the 
corresponding momentum integral. 

It remains to solve Eq. (||). In the case of non-exceptional 4-dimensional kinematics, rank(S) — 
rnin(N, 6) holds. For N < 6 one has k = -1/2 J2k=i Ski 1 - In the 

case iV = 5, contains a factor 
(N — 1 — n) which is 0(e) . As is well known, pentagon integrals are just a sum of box integrals up to a 
remainder which drops out in phenomenological applications. In the case N > 6 and 4-dimensional 
kinematics, Eq. (^) does not have a unique solution. To clarify this point it is useful to rewrite 
Eq. (ph. Using momentum conservation, rjy = 0, Eq. rtq) leads to the following equations 



JV-l 



JV-l 



JV 



^ Gkibi = v k B N , 2_j Vlbl = 1 ' B N = y^ k 



(11) 
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Herein G is now the Gram-matrix of the vectors n=i...jv-i- In the case of four dimensional kinematics 
it is at most of rank 4 for all N > 5. Generically any four vectors out of ri, . . . , rjv-i are linearly 
independent. Such a configuration is called non-exceptional in the following. Choosing four linearly 
independent vectors Eft. 1 4 as a basis of the physical Minkowski space, one can define a coefficient 
matrix R and a Gram matrix G made out of these basis vectors 



^ R ij E i i G jk =2Ej- E k (f 2) 

i=i 

The Gram matrix G is expressible as G = R T G R. 

Based on (|l2|) one can construct the most general solution for the case det(G) = by means of 
pseudo-inverse matrices. This concept will also be useful for the tensor integrals. 

A pseudo-inverse H to G is defined by the conditions H G H = H and GHG = G. Given a 
pseudo-inverse matrix H to G, the following statement can be proven. A solution to the linear 
equation G ■ x — y exists if and only if y = G H ■ y. Then the most general solution can be written as 
x = H -y+ (ljv-i — H G) -u, where the last term — with u € R N ~ X arbitrary — is parametrizing the 
solutions of the homogeneous equation. As G is symmetric, its pseudo-inverse is uniquely defined 
M . The pseudo-inverse to G is given by 

H = R T (RR T )- 1 G- 1 {RR T y 1 R (13) 

where R is defined in ( |l2| ) . For N = 5 and non-exceptional external momenta the vectors r% , . . . , ri 
are the natural basis. Then R = I4, G = G, H = G _1 . Note that the concept of the pseudo-inverse 
always allows the inversion of linear equations containing the Gram matrix, even if the external 
momenta are exceptional^. 

If the Gram matrix is not of maximal rank, the condition v — GH ■ v for the existence of a 
solution of the inhomogeneous equation in (|ll| ) is never fulfilled. This can be seen by diagonalizing 
the symmetric matrix GH by an orthogonal transformation O, G H — O t DO, where D is a 
diagonal matrix with elements (1, 1,1,1,0,..., 0). Now it is evident that in general O ■ v ^ DO ■ v. 
Thus a solution to Eq. ( ]li"| ) exists only for the case Bn = 0. The solution of the homogeneous 
equation spans a (N — 5)-dimensional space which is just the kernel of the Gram matrix. It can be 
parametrized by (N — 5) vectors [/C 1 '---'^ -5 ). Defining 

K = l N _ x - HG = ljv-i - R T (RR T )- 1 R , (14) 

one has K ■ v 6 fcer(G). Now one can choose U^ N ^ 5 ^ = K ■ v/[v ■ K ■ v) parallel to v and the 
others orthogonal, v ■ J7( fc=1 >— i^ -6 ) = 0. A general vector in ker(G) is then parametrized by 
U = J2k=i PkU (k) + aU {N -^ and the general solution is given by 

bt _ + E^M>f> , = 1 ,... >A ,_ 1 

v ■ K ■ v 

N 

B N = ^6 fe = (15) 

k=l 

where a — 1 is imposed by v ■ b = 1 . Thus for N = 6 the solution to ( |TT| ) is unique and therefore equal 
to the one defined by the inverse of the matrix S. Eq. ( |l5| ) proves constructively that for all N > 6, 
I N can be expressed in terms of (JV — l)-point functions without the higher dimensional remainder 

2 Clearly one could also solve Eq. ([J) by means of the pseudo-inverse to S. With Os defining the orthogonal 
transformation which diagonalizes S and the non-zero eigenvalues Ai , . . . , \q of S, the pseudo-inverse to S is given by 
Og • D'g lv ■ Og, where D l g v = diag(l/Xi , . . . , l/Ag, 0, . . . , 0). Its computation is algebraically more involved than the 
computation of H though. 
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term. This is a consequence of the linear dependence of propagators if the external momenta are 
4-dimensional. For the special case JV = 7 this has already been demonstrated in B. Obviously one 
can choose the (3's to eliminate (JV — 6) 6's from the set {bi, . . . , b^-i}- Doing so one observes that 
I N can be expressed by only 6 (JV — l)-point graphs for arbitrary JV > 6. Here we make contact to a 
result in |j, 18 , where a similar relation has been derived for the IR finite case in integer dimensions 
by using 4-dimensional Schouten identities. 

For JV < 6 one has det(5') ^ and the following relation (see also 0) holds 

1=1 l,k=l y ' 

which shows that the vanishing of the finite remainder terms (|l^) is related to the vanishing of the 
Gram determinant. 

In summary, with the definition of reduced graphs 

jn _ f ] rim=l iL jn _ f , IIm=l (gjV ~ %n } (m 

1 N-pM,-,l P - / aK n JV 2 ' 1 N-p,N-h,...,N-l p - / aK , K 1 1 ) 

the scalar reduction formula for a regular matrix S is 



1 N 



N 



\ E S u r N-U + (N-n-l) ^f^y I n N +2 , det(5)^0. (18) 



In the case of a singular matrix S, as outlined above, we can always achieve a representation without 
higher-dimensional integrals. 

N-l N-6 

In = -^ 7 J^T,(( K -^+T,^ U i k) ) I N^N-i , det(5) = 0. (19) 

1 = 1 k=l 

Remember that the (N — 6) /3's are free parameters and that the jy( fc=1 >---,iV-6) are lying j n the 
kernel of G and are orthogonal to v. The case N = 6 is special in the sense that it is of the form 
(18) with vanishing higher-dimensional term. 

By applying the above formulas iteratively, any iV-point function can be reduced to linear combina- 
tions of box integrals plus irrelevant 0(e) terms in a constructive way. 



2.2 Application 

We will now discuss the formulas for the cases up to N = 6 in more detail. In these cases one has 
det(5) ^ and relation ( ^6|) holds. Note that the vanishing of the Gram determinant for N < 6 
typically occurs at a border of the respective phase space. In the case N = 3 (4, 5) and det(G) = 
the scalar integrals are just sums of 2-point, (3-point,4-point) functions respectively, an observation 
made some time ago by Stuart |l|| and reflected by Eq. (|l8|) . This fact can be used as a guideline to 
define nontrivial groupings of (JV — 1)- and JV-point functions which is helpful for numerical purposes 
[17]]. The other dangerous determinant is det(S') which vanishes if JV > 7. We note that in this case 
our formulas are mathematically well-defined and that for the generation of explicit expressions one 



just has to apply formula (19). Due to the freedom to choose the parameters /?, the representation 
of the JV-point function in terms of reduced integrals is not unique as it is the case for JV < 6. This 
is a reflection of the fact that the respective reduced integrals are not linearly independent. 

As the reduction formulas were derived under the assumption JV > 4, we discuss now the special 
cases JV ~ 1,2,3 first. We will represent the kinematic information in terms of the matrix S. The 
relation to the Gram matrix is defined in Eq. (|^). 



■5 



Cases N=l,2: 

As massless tadpole integrals are zero in dimensional regularization, /" = 0. The case N = 2 is 
trivial in the sense that reduced graphs are tadpoles. The formula simply gives a relation between 
n- and (n + 2)-dimensional two-point functions. If p\ is the external momentum with p\ = m\ , one 
obtains 

/« = 2^-—^ 7 2 "+ 2 (20) 
For vanishing Gram determinant, i.e. m\ = 0, = in dimensional regularization. 
Case N=3: 

The matrix S generally looks like (pj, - 

1 / m\ m\ \ 

n\ ml (21) 
n\ m\ J 

Note that in the cases of one or two external lines on-shell, S is not of maximal rank. If for example 
m\ = m\ = 0, or m\ = 0, m 2 , 7^ m§, there exists no solution to (||) and no reduction is possible. If 
all legs arc off-shell one finds a solution which relates an (n + 2)-dimensional off-shell triangle to an 
n-dimensional one. 

I3 = -- T S^I?, + det{G \ (2 - n)/"+ 2 (22) 
3 2 j-± lk 2 ' 1 8dct(S*) V ' 3 V ; 

We will see later that the (n + 2)-dimensional off-shell triangle appears in the tensor reduction of 
rank > 3 tensor 4-point functions. 

As any A-point integral in (4 — 2e) dimensions can be reduced to triangles plus a finite remainder, 
we see that they are the atoms of any scalar reduction formula. Moreover, at the one-loop level they 
define an IR counterterm structure to any A^-point function)^]. 

The explicit expressions for the three types of triangle graphs may be found for example in |l7| . 

Case N=4: 

With (pi +P2) 2 = s, (p 2 +P3) 2 = t, one has 
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m\ \ 
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ml 
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m\ 
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ml 


/ 



S = -- "r "o "., (23) 

S is of maximal rank for all on-shell cases as long as s and t are non-zero. The reduction formula 



(18) applies and relates a 4-point function in n dimensions with triangles and (n + 2)-dimensional 



4-point functions. 



Lk—1 



For n = 4 — 2e the triangles carry all the infrared poles whereas the (6 — 2e) dimensional part is 
a finite remainder. It can be calculated directly by setting e = 0. Explicit expressions for the box 
integrals can be found in |l7j . 

3 This is actually to be expected from looking at the reduced diagrams (obtained by shrinking all finite propagators, 
if the loop momentum becomes soft/collinear) of the one- loop iV-point function corresponding to solutions of the 
Landau equations. The respective reduced diagrams are just the reduced diagrams of triangle graphs. 
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Case N=5: 



With (pj +Pj+i) 2 = Sj,j+i (j mod 5) one has 



S=-~ 
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For = 5 the higher dimensional integrals come with a prefactor (4 — n) 
phcnomcnological applications. One only has to know box integrals. 



(25) 



and thus drop out in 



T n — - 



1 
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Lk=l 



C-l rn 
D lk 2 4,L 



0(e) 



(26) 



Explicit expressions for the 5-point function with zero or one massive legs may be found in 
Other cases are easily constructed by using known representations of off-shell box- integrals. 



Case N=6: 

Now the Gram matrix is not invertible anymore for 4-dimensional kinematics whereas an inverse for 



S still exists. With (pj + pj. 
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(j mod 6) one has 



(27) 



Before going on, it is instructive to clarify the dependences between the Mandelstam variables which 
define the matrix S. Since they are directly related to the cuts of the Appoint graph, the following 
counting holds for TV > 4, where M is the number of off-shell legs, C the number of cuts, D the 
number of independent Lorentz invariants built out of the vectors ri,...,rjv, and Z = C — D the 
number of constraints: 



C = N(N -3)/2 + M , D = 3N- 10 + M , Z = (N — 4)(JV — 5)/2 



(28) 



This makes manifest that for N > 6 one encounters subtleties due to the appearing constraints. Here, 



the (nonlinear) constraint det(G) = shows up in the fact that X^ti h 
The reduction formula for N = 6 reads 



1 v^JV o-l _ r, 

2 Z^fc,2 = l °fci — u - 



1 6 

r« — 9- 1 r 



— 1 

5.i 



(29) 



l.k=l 



If all six external legs are on-shell, the occurring pentagon integrals Ig l have one external leg off-shell. 
The explicit expression for the on-shell 6-point function is given in Appendix [A|. 



3 Reduction of tensor-integrals 

In this section we show that any rank L, Appoint tensor integral can, by recursion and scalar 
reduction, be expressed in terms of scalar box, triangle and two-point integrals. We prove that for 
N > 5 and non-exceptional external momenta higher dimensional integrals drop out. 
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The rank L, A-point tensor integral is given by 



k^ ... k^ 



(30) 



We omit the superscript n indicating the dimension in the tensor integrals to simplify the notation. 

If the dimension is different from (4 — 2e), it will always be written out explicitly. 

After introducing Feynman parameters and making a shift of the loop momentum as in the last 

section, the odd powers of the loop momentum in the numerator can be dropped. Using the standard 

integral 

{/il.../i2m} f 1' 

2. 



dn- 



k^ . . . k^ 2 " 
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'(k 2 -M 2 ) N 
one finds the following formula: 



9{m) 
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1 N 
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m=0 



Jl, 



N-l 
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• jL-2m = l 



9(m) r j 



T(N) 



{Hi...Hl} 



M 



-N+(n+2m)/2 



(31) 



,3L-2m) (32) 



[L/2] is the nearest integer less or equal to L/2 and 



9'(m) r h 



3L-2n 



{Hi—Vl} 



stands for the sum 



over all different combinations of L Lorentz indices distributed to m metric tensors and (L — 2m) 
r- vectors. These are ) IlfeLi(2fc — 1) terms. A dot appearing as an index at objects inside a 
square bracket stands for one index out of the set specified in curly brackets at the outside of the 
square bracket. iW or Xf m \ denotes the product of m terms of X with adequate indices^]. For 

{A 1 1A 1 2A 1 3/J4A'5} 

with X a symmetric tensor of rank 2 and Y a 



example, a tensor object like 
vector means 

{/JlP2M3A 1 4A t 5} 



x vi^ x ^^)Y^ + 5 permutations. 



■ ■ ■ ,3m) are scalar integrals with nontrivial numerators in D dimensions, defined by 



In til, ■ ■ -Jm) 



(-1) N T{N - D/2) I d N z5{\ 



N 

E 

i=i 



°3l 



JM 



{z-S-z) N - D / 2 



(33) 



Recursion relations for this kind of integrals were given for M < 4 in [[i~7| . We derive such relations 
for general M in Appendix |b| 

From now on we use standard summation conventions for the indices. Eq. (|32| ) implies that tensor 
integrals in momentum space are linear combinations of scalar integrals in different dimensions. We 
note that formula ([52|) is equivalent to a more general formula given in jllj] for the case of integer 
powers of the propagators. 

We want to prove now that for non-exceptional momenta and N > 5 the higher dimensional 
scalar integrals drop out. To this end we solve Eq. ( |32| ) for I N (ji, ■ ■ ■ , j'l) by contracting it with 
2 L rf^ . . . r[^ L . For the inversion of the Gram matrices we use its (pseudo) inverse as defined above. 
With 2ri ■ k = q 2 N — qf + vi one finds 

/vU: = EK~ fc) 



k=0 



H, 



■ H, 



N-k,N-l u ...,N-l k 



[L/2] 

E 

rn — l 



(-iy 



(m) r n+2m 



mr'i 



N 



((£ — 2m) indices) 



{ji,— Jl} 



(34) 



4 The "power" (m) may appear as a lower index for convenience of notation if it would interfere with the dots 
standing for upper indices. Hence we write X (m) instead of (X'-)C m ). 



where w — H ■ v and H is defined in Eq. (|l3|). On the right-hand side of Eq. (|34|), differences 
of reduced integrals appear as defined in Eq. (|l7]). In principle one has to add to j34| ) also the 
solutions stemming from the homogeneous equation which are present if G is not of maximal rank. 
But as they contain the matrix K = In-i — HG, they vanish after contraction with the r's, 
i.e. ■ (1jv-i - H G) = 0. If G is invcrtible, H = G~ l . The bracket of the fc-th term of the 
first sum stands for (^) terms whereas the bracket of the m-th term of the second sum stands for 
Gm) n&Li (2fc — 1) terms. Inserting ( |34| ) into (|32|), we find that the tensor integral decays into a 
part containing (4 — 2e)-dimensional objects, K^" , and a part built out of higher-dimensional 
objects, J^- Mi . 



N 



J/J.1.../J.L 



{lll,...,llL} 



N ' J N 

L 

^ [ W (L-fc)^h ■ ■ -^h ^N-k,N-h,...,N-l,, 

[L/2] N-l 

E(- 1 ) m E [{^m- {m) -u (m) y n 

m=l Jli...,3i-2m=l 



(35) 
(36) 



{Ml— M-E,} 



(37) 



The Lorentz indices are carried by the objects 





= r"- 


H -r 




= (r^ 


■H)t 


w 


= r^ 1 • 


w = 



KP-v 



(38) 



We recall that w — H ■ v and vi = Gu/2. 

If the r's span 4-dimensional Minkowski space, which is generically the case if N > 5, TL^ is 
just proportional to the metric tensor g% in 4 dimensions: 



N-l 



4 

E 



E 



! ^k^lk - 



(39) 



Thus the coefficients of the the higher dimensional integrals in (^) are proportional to — g± 
which is of order e. As by power counting the higher dimensional integrals are finite objects if the 
external momenta are 4-dimensional, the whole contribution is 0(e). This proves the cancellation of 
higher dimensional integrals in tensor reductions for arbitrary N > 5 and non-exceptional external 
momenta, and we obtain 



n 

k=0 



(L-k) 



I 



N-k,N-l u ...,N-l k 



0(e) 



N > 5 



(40) 



For N < 5, Eq. ( p9| ) is not valid since the external momenta cannot span Minkowski space anymore. 
Thus one has to calculate the terms J^'"^ L for N < 5. The explicit expressions are given in 
Appendix [E)|. 

Now we want to rewrite Eq. ( |35| ) as a recursion formula for arbitrary tensor integrals. To this 
effect we express the contracted tensor integrals I N _ k jy-jj N-i k s • > 0) a s (N — l)-point tensor 
integrals which are maximally of rank (L— 1). By using (k— 1) times the relation qjj — qf = 2r;-fc — 
one gets 



rn 

1 N-k,N-L 



,N-l k 



k-1 

E(-d" 

3=0 



j-1 



V 



v: ' r., 



{h...i k -i} 



N-l, N-l 



-u,k>l (41) 
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Insertion of expression (Elf) into Eq. (pq) leads to the recursion formula 



TfHl...llL _ 

N - 



WK 



r {L-l dots} 



{fll...HL} 



2< 



L-l) 



LI 



l N-l,N-l 



(42) 

The derivation of (f42| ) is given in Appendix ^|. 

Since for A > 5 the terms J^ 1 ' contribute only to order e, K^"'^ L in Eq. ( |42] ) can be replaced 
by I&'"^ L to obtain the recursion formula for A > 5. For A < 5, one cannot drop the J^"'^ L 
terms, such that the general recursion formula reads 

1 



1 N 



TfJ.l...fJ.L 

N 



W{I N - J N ) {L ~ ldots} 

{Ml— Mi} 



■ K, 



L N-l,N-l 



(43) 



As a consequence, higher dimensional integrals I^t 2m (m — 1,2) will appear in the reduction of 
tensor Appoint integrals only with A' < 4. All the higher dimensional integrals can be re-mapped 
to n-dimensional integrals with the scalar reduction formulas (|20|), ( p2| ) and (|2^). 
In the case of exceptional kinematics and A > 5, higher dimensional (A > 5)-point functions can 
be present. 

As some of the reduced integrals on the right-hand side of ( f43| ) do not contain a trivial propagator, 
they are not in the standard form for applying the reduction formula again. Therefore one has to 
perform a shift k — > k + ri in the tensor integral, leading to 



1 N-l,N-l 



(R) 



1 N-1 



(R[N]) 



l N-l 



k=0 



•(fe) r {P-k dots} 



(R[i\) 

{Hi—Lip} 



(44) 



The argument vectors are 

R 

R [k] 



R[i] - 



(n, ■■■ ,r N ) 
(ri,...,f fc ,...,r N ) 

(r ;+ i - n,n + 2 -n,. 



-,r N - 2+ i -n,o) 



(45) 



where f means that the respective vector has to be left out. The vector indices are understood to 
be taken cyclically symmetric with periodicity A, i.e. tn+i — n for I S {1, . . . , A}. As we assume 
rjv = 0, the integrals in the last line of (|44|) have again at least one trivial propagator and are suited 
for a further reduction step. 



A full tensor reduction of a rank L tensor integral is obtained by employing Eq. (42) for N > 5 
resp. Eq. ( |4"3| ) for A < 5 to do the first reduction step. Then the integrals iN—itf—i have to be 
shifted by using Eq. ( pl4| ) to end up with expressions of rank < L — 1 to which the same reduction 
procedure can be applied again. Explicit tensor reduction formulas for N — 2, ... ,6 are given in 
Appendix [j5|. 

For A > 5 one can express by recursion any rank L, A-point tensor integral in terms of (4 — 2e)- 
dimensional scalar TV-point and tensor box integrals, which in turn can be reduced further down to 
scalar box, triangle and two-point integrals, taking also into account the scalar reduction formulas 
given in the previous section. This defines an algorithm which can be easily programmed. 
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4 Conclusion 



We have presented reduction formulas for massless Appoint scalar and tensor functions. 

We have shown in a constructive manner how scalar Appoint functions can be represented as 
linear combinations of (4 — 2e)-dimensional scalar (N — l)-point functions for arbitrary N > 5. 
In particular, we pointed out how to treat the linear equations which determine the reduction 
coefficients, bi, for N > 6 in terms of the pseudo- inverse of the singular Gram matrix. In this way 
all mathematical operations are valid for 4-dimensional external kinematics throughout the whole 
reduction procedure. 

We applied the reduction formalism to scalar integrals up to N — 6 explicitly and gave an expression 
for the on-shell 6-point function. 

For Appoint tensor integrals we have formulated a reduction scheme for arbitrary N. We have 
proven that higher-dimensional integrals always vanish for N > 5 in the case of non-exceptional 
kinematics. We derived a recursion formula for the remaining n-dimensional part. For the derivation 
we used methods a la Passarino-Veltman, such as contracting tensor integrals with external vectors 
and inverting Gram matrices. In the case of a singular Gram matrix the inversion has to be done 
with its pseudo-inverse. By iteration any tensor integral can be expressed as a linear combination 
of scalar integrals. Higher dimensional scalar integrals appear only in the reduction of N < 4 tensor 
integrals. These higher dimensional integrals are expressible in terms of (4— 2e)-dimensional integrals 
as has been shown explicitly in the discussion of the scalar reduction formulas. To make the general 
formalism more user-friendly we gave explicit expressions up to N = 6 in an appendix. 

We also derived a reduction formula in Feynman parameter space in an appendix. The formula 
generalizes results in the literature avoiding a projective transformation. As we gave all the formulas 
to translate momentum space expressions into Feynman parameter space expressions, the equivalence 
between the two approaches is manifest. In applications this will allow us to employ methods similar 
to the ones in |l7j to get numerically stable expressions, a problem we did not address in this paper. 

We conclude that the computation of IR divergent one-loop integrals for arbitrary numbers of 
legs can be mastered with the reduction formulas presented here. The iterative structure makes it 
easy to implement the formalism in algebraic computer programs. The conceptual problems for the 
construction of multi-parton one-loop amplitudes are thus solved. 

The generalization of our method to include also massive particles is postponed to a future 
publication. 
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A The scalar 6-point function 

Here we present an explicit formula for the scalar hexagon function with all external legs on-shell. 
The reduction formula (p9h allows to write it as a sum of six pentagon integrals with one external 
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leg off-shell. 



^6 l ( S 12i S 23, S34, S45, S56, Sqi, S123, S234, S345) 



1 6 

- y 

2 ^ 

fe,;=i 



5 , fc; 1 ^5,lmass(s;+2,i+3j 5/4-3^+4, S; + 5^f + i, 5;j + ij + 2; 5;j + i) 



(46) 



The function I5.ima.ss has been given in We rederived the formula. Although we agree with 
it for Euclidean kinematics we note that the expression given there is not well suited for analytic 
continuation 5 . The other ingredient, the matrix 5, is given explicitly by Eq. ( [27] ) after setting all 
mj = 0. The indices in ( ]46| ) have to be taken modulo 6, and due to momentum conservation one 
has S456 = S123, S612 = S345, S561 = S234. The nine Lorentz invariants are not independent since 



they fulfill the constraint det(G) = 0. We find 



Ie = 

Ai = 

B\ = 

Ci - 



6 

E 

k=l 



i-^12)" e 



B k ) + C k + D k 



(47) 



S61S12S23S234 

»X2 



(Sl2 - Sl23) 2 [(-Sl2)" 



-S123J 



(Sl2 - S 345 ) 2 [(-Si2) e - (-S345)" 



S12S23S34S123E1 
I ^ S45 — S345 

S34S45E1 

-s 12 f 



S34S45S56 

(Sl23 



Dl = -77 



+63 

63 

+61 

62 

+b. 



S12S23S34S123-E1 
(5345 ~ S12) 2 

S56S6lSl2S345^1 
1 , S12 



Li 2 1 - 



S45 — S123 
S45S56-S1 

5l2 
S123 
5l2 
S345 
Sl2 - 



Li a 1 - 



S56S6lSl2S345^1 
S123S345 



- Li 2 1 - 



Li 2 1 - 



S123 



-Li 2 

S345 



1 - 



S12S45 
S123 
S12 
S345 
S12 



S34S45S56 3 



S61S12S23 
1 

S56S61S12 
1 



+ 



S12S23-E1 S61S12-E1 

561 — S234 Hi ~ s 345 



S56S6I-E3 
1 



556 s 6l5l2 S56S61S345 
/ Sl2 — S345 S12 — S123 



\S5eS345E1 
1 

h 

S45S56S6I 

1 

S12S23S34 
512 ~ 5123 
5345123^1 



S34Sl23^1 
g56 ~ 5 123 
S45S56-E'2 
, S23 — S123 



S6IS12E3 
1 

S61S12S234 



1 2 / s 12 

log 

S 6 i 



S234 S 56 — SX23 



556 



V Sl2 5234^2 
S234~ 



1 



S12S23-E2 

-b 4 



S56SQlE2 

S23 — 5234 
S23S34E2 
S12 - 



log 2 



S34S123E2 



5123 
S234 



log 2 



£12 

S56 



+ 



s 345 



S61S12S23 seiSuEi 
S56 — S123 , S23 — S123 



1 2 / s 12 

log 

S123 



1 



1 



5615125234 556561512 



,5343123-^2 S12S23-B2 

561 — S234 S34 — S234 \ 
S56561-B3 S12S234E3) 

5 The problematic term (up to a trivial typo) is the bracket (1 — (m|s23/s4s/^5l) — e )- Rederiving the formula we 
actually get (1 — (mg/s45) _E (s23/s5l) — e ). Only in the latter form the analytic continuation is defined by simply 
adding an infinitesimal imaginary part to the Mandelstam variables. Care has to be taken as well in the analytic 
continuation of Lt2(l — "15^23/545/551), cf. the remarks at the end of this Appendix. 



12 



bi 



1 



1 



S23 — S234 — S234 



S12S23S234 
/ S12 - S345 

V S56S345-K1 



S12S23S34 
1 



b 4 



S23S34E2 S12S234E2 
S12 - S123 



log' 



(-) 

\ s 234/ 



S6IS12S23 



+ 



with the abbreviations 



S12S23-E-1 
s 61 — s 345 



/ S34 - S345 _l_ 

\ s 56 s 345-£'3 S6IS12E3 



E 1 


= E 4 


— S123S345 


- S12S45 


E2 


= E 5 


= S234S123 


- S23S56 


E 3 


= E e 


= s 345 s 234 


- S34S6I 



log' 



(-) 

\ S345 / 



and 



r r = 



r(i + e)r 2 (i-e) 



T(l - 26) 

The functions Ak,Bk,Ck, -Dfc for k > 1 are obtained by cyclic permutations of the indices. The other 
necessary ingredients are 

7 / 2 2 1 2 

01 — ( — 5345555523 + S12S234S45S34 — S56 S 23S345S45S34 — S45S234S345Sl23 s 34 

+S34S45S6IS123 — S34Sl23 s 12S234S45 ~ S34S123S56S23S345 + 2S34S45S23S12S56 



+ s 34Si23S345S234 



2 2 
S 34 S 61 S 123 



«34S345 S 56 S 123S61 



1 2 1 2 

— S45 s 234 s 123 s 56 s 345 + s 45 s 234 s 12 s 56 + s 56 s 23 s 345 s 45 



s 345 s 56 s 123 s 234 

— •S56 s 34 s 45 s 61 s 123 



+2S34S12S56S61S45 — S12S56S45S234S345 ^ S12S234S45 ~~ S 123 S 234S345 

+ S123S56S23S345 + 2S45S345S234S123S12 ~ 2S34S45S12S56 + 2S45S34S123S56S345 

— S4 5 S34S 6 lSl23 s 12 — S45S56S23S345S12 + S34Sl23 s 6l5345)/-F 

where the other o^'s are again obtained by cyclic permutation. Finally 



(48) 



F = 



— s 



:S45- , 56°23 



4s23S34Si 2 S5 6 S 6 iS45 - 2S23S34S345S56S123S6I ~ 2S23S12S56S45S234S345 



+ 2S23S 3 4 5 S56S123S234 — s 123 s 61 s 34 — 2s34Si25i23523456l545 + 2s34S 123 S6lS345S234 



S A rS 



45*234*12 



s 345 s 123 s 234 + 2si2Sl23S234 s 45S345 



64 det(S) 



(49) 



Note that in the e-dependent part of the formula no 6's appear. We did not succeed in finding a 
more compact form for the part which does not depend on e without spoiling this nice feature. 

The expression for Ig in the form as given above is strictly only correct in the Euclidean region 
where all Mandelstam variables are negative. For most of the terms, the analytic continuation to 
positive values is defined by simply using the replacement s — > s + iS, where s stands here for any of 
the Sij, Sijk- No cut will be hit by the logarithms, the dilogarithms with a single ratio of Mandelstam 
variables and the exponentials (— s — i5)~ e . Concerning the dilogarithms of a product of ratios, more 
care has to be taken. To avoid the crossing of a cut one has to make the replacement 



Li 2 1 - 



5152 
5354 



Li 2 1 



si + iS s 2 + iS 



+ V 



S3 + iS S4 + iS 
si + iS s 2 + iS 



v S3 +iS' S4 + iS 
\og{xy) - \og{x) - log(y) 



log 1 



s\ + iS s 2 + iS 
S3 + iS S4 + iS 



(50) 
(51) 
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B Recursion relation for integrals with Feynman parameters 
in the numerator 



Here we want to comment on the reduction of Feynman parameter integrals with nontrivial numer- 
ators as defined in (|33|). In recursion relations for integrals with up to 4 Feynman parameters 
are given. The derivation is based on the approach of J3j which uses a projective transformation 
|2C|| . We present an independent derivation and generalize the formulas to arbitrary N and numbers 
of Feynman parameters in the numerator. 

Consider the following identity (j — 1, . . . , N — 1) 



7 r) ( N N \ 

j dNz ^- [ n - e **k ■■■ziAz-s- Z )- N+n/2+i j 



(52) 



Setting tn = 0, one obtains z ■ S ■ z = Y2kl=i z iGikZk/2 — Sltl] 1 v k?k- The first step consists in 
eliminating the (^-function by integrating out z^. If now the derivative acts on 6- functions, terms 
with (5-function insertions are produced. The reduced integrals I^-i j (defined in Eq. ( j53|) below), 
obtained by pinching the j'-th propagator line in an ./V-point graph, correspond to these (5-function 
insertions. If the derivative acts on the Feynman parameters, the monome in the numerator is 
reduced by one degree, and if it acts on the z ■ S ■ z term it formally decreases the dimension by 
two and increases the numerator by (G • z)j — Vj. After these operations one reintroduces the z/y- 
integration with a delta- function insertion. Using the following generalized definition for pinched 
scalar integrals with nontrivial numerators 

I%_ hj (h, ...,l p ) = (-l) N -iT(N 1 - n/2) r d N zd(l - £ z t ) \ ' (53) 

Jo l=l {z ■ o ■ z) i 



we find 



N-l 



E GjIoInQo) ■ ■ ■ , lp) = ^ fi]i k lN +2 (h, ■ ■ ■ , Ik, ■ ■ ■ i P ) + In-i,n(Ii, ■ ■ ■ Jp) 



k=l 



-I^_ hj {h,...,l p )+v j lUh,---,l P ) (54) 



Herein Ik means that the respective index does not appear. 
Now we want to sketch the proof for a second equation 



N-l 



I%_ ltN (h, ...,l p ) = (N-n-p- l)I^ +2 (lu ■ ■ ■ , l P ) + J2 v i In(Io, ...,l P ) (55) 

The proof is done by induction to p. In order to show that the induction start (p = 0) holds true, 
one directly calculates 

1%-iA 1 ) = / dK T^~ 2 =(N-n- l)i£ +2 (l) + E v lo I%(l ) (56) 



Z = l 



For the induction step, one assumes that (|5J]) is fulfilled for a given p. Viewing the vj as independent 
variables, one differentiates the formula with respect to vi p+1 and finds the formula for p + 1 and 
n — 2. As the dimension is arbitrary we can replace n ~ 2 by n in the expression. This proves the 
validity of equation (J55|) for all p. 
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Combining (|54|) and (|5q) one finds (j = 1, . . . , N) 



N 



2 E ^Ho^nO'O, ■■■Jp)—/ j fiji k lN +2 (h, ■ ■ ■ , 4, ■ ■ ■ , Zp) 



(57) 



io=l 



+ {N-n-p- l)I^(h,. . . , l P ) - I^ ld (h, ...,l P ) 

In the case iV < 6, S is invertible and the final reduction formula reads 

p 

2 



1 

I 7 ^(Iq, . . . , l p ) — — ^ S lo ] k I r ^ 2 (li, ...,lk,...,l p ) 

k=l 

N N 

+o E sifo 1 ^ - » - p - w^ih, . . . , g - 5 E ^o^-i,i(^' • • • > w 



(58) 



This is the generalization of the formulas given in |i7| for the case of four Feynman parameters in 
the numerator (p = 3). (In their conventions S7L = riik&iak/N n and '^2 i S~ k 1 = 7fc/-/V„, which is 
identical to our S 1-1 up to a trivial relabeling of indices). It means that all Feynman parameter 
integrals with numerators can be reduced to ordinary scalar integrals by iteration. 

In the case N > 6, 5 is not invertible. The linear dependence of the row (resp. column) vectors 
of 5* makes it difficult to generalize the Feynman parameter space based techniques to arbitrary 
N H [l^] . In our approach one can use Eq. (|57]) directly in this case, where the inversion of the 
Gram matrix should be done with its pseudo-inverse as explained in the main text. As already noted 
earlier, it is also possible to work with the pseudo-inverse to S. In any way, as we have proven, it is 
the (pseudo) inverse to the Gram matrix which induces the cancellation of the higher dimensional 
terms by virtue of Eq. ([39]) . 

Finally, with the formulas given above and in the main text, it is possible to translate expressions 
in Feynman parameter space into expressions in momentum space and vice versa for arbitrary N. 



C Derivation of the recursion formula for tensor reduction 

In this appendix we give the derivation of the recursion relation (E2I) : 



7^/il...A«L _ 

N 



■ jr{L — l dots} 



WK 



N 



{Mi-Mi} 2( L ~ 1 ) 



{aii... pt} 



N — l,N — l 



(59) 



To Eq. © 



"TV 



{Ml— Mt} 



pi 

1 N-k,N-l 1 ,...,N-l k 



we apply Eq. ( [ill ) 
/ 



N-k,N-h,...,N-l k 



3=0 



3- 



(k-j-l) 

V} r. Ul . . . r. v 



{h-h-i} 



I 



N-l,N-l k 



and get, by taking into account combinatorial factors from contracting the tensor brackets 

l , fc-l 



Tf/J,l...fJ,L 



fe=l j=0 ^ 



W(l~}-\)K{K Vi ■■■'H V] 



N—l.N—l 



(60) 
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Using the fact that 

L k-l 



L L-l L-2 L 

EE a ^ = EE 6 tt - k ~ l ) ak i = a L,L-i + Y ak i 

fe=l j=Q k=l j=0 j=0 fc=i+l 



(61) 



we can separate terms with and without W-vectors. 

L L-2 



... +EE^'^ i; - 1 ) 



2*(_l)*-i-i / L-j-l 



fc=l j=0 



fc - i - 1 



(62) 



{Mi—Mi} 



7 



N-l.N-1 



,L-1 r 



L! 



?d ■ ■■U„ L _ 1 K l 



{Mi—Mi} 



1 N — l.N — l 



Now one has to rearrange the terms containing W-vectors. 

First we write 1/fc in ( |62| ) as 1/fc = (Z, — k)/{Lk) + 1/L. The 1/L part vanishes since the sum over 
fc just represents (1 — fp-i— l = 0. The (L — fc)/(Lfc) part makes the upper bound of the fc-sum to 
be L — 1. From the tensor bracket a W can be factored in a symmetric way: 



{Mi— Mi} 



W Ml 



L-j-l 



W( i _ J _ 2 )W„ 1 . ..K Vj fC\ 



■ Ti'.K) 



{hi...hl-i} 



Hence Eq. (B2l) is equivalent to 



fc'Ml— Mi 

N 



-)L-l r 



LI 



Hi 



{Mi— Mi} 



l N-l,N-l 



w m 1 ^ ^— ; 



L-l fc-l 



2 j (-i)k-J-i / L-j-2 



k=i j=a 

|{M2— Mi} ,., ... 



fc - j - 1 



[L — 1) permutations 
Comparing with Eq. (^30|) we see that the right-hand side of Eq. (^33|) is just 



- (L — 1) permutations = — 

L/ 

Combining (^) and ([54]) we obtain the recursion relation ( j59|) 



^.^{L-ldots} 



{mi-.-M-c} 



(63) 



(64) 



D Explicit reduction for N=2,3,4,5,6 

We now give explicit formulas for tensor integrals up to N = 6, and rank L < N. They have 
been derived by applying Eqs. (f43|), ( pfjj ) and (|37|). The reduction stops because tadpole integrals 
of massless propagators are zero, I\,I X = 0. This shows that any tensor integral with an arbitrary 
number of legs can be expressed in terms of scalar integrals only. For non-exceptional kinematics 
higher dimensional integrals (m = 1,2) occur only with N' < 4 in the reduction of arbitrary 

iV-point tensor integrals, as has been explained above. All the higher dimensional integrals can 
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be re- mapped to n-dimensional integrals with the scalar reduction formulas (J20|), ( |22| ) and (E 
Explicitly, the expressions for the higher dimensional integrals are given by 



TMlA»2 

,7 JV=2,3,4 

TAtlA»2A»3 

J Af=3,4 



J, , . . . 1 {M1M2M3} 

(-l)[( 5 /2-W) Wj 

(5/2 - W) (2) 
1 



(65) 



rn+2 
J 7V 



(g/2-ny^ 



{Atl^2Al3} 



n+2 
'JV-1JV-Z 



T-n+4 



( 5 /2 - W)" W, 



(2) 



{^1/^2^3/^4} 



(5/2 - H)" W/Q 



{A t lM2M3A t 4} 



-n+2 
'3,4-Z 



(g/2-nyn^ 



T n+2 
J 4 

{/il "2M3/J4} 



f,n+2 
'3,4-Z 



Note that rank one, (n + 2)-dimensional three-point functions appear. As the reduction rules do 
not depend on the space-time dimension all the formulas given above are also valid for higher 
dimensional tensor integrals. The tensor structure is carried by the well-defined 4-dimensional 
objects W , fCf , given in Eq. (|38|). Note that metric tensors occur only in the combination 
(g/2 — H) as coefficients of the higher dimensional integrals 7JJ- m . 

For N > 5, the equal signs are only valid up to 0(e) since the Jjt J 1 " 41L terms have been dropped. 
The shift of the integrals InZtjv-i(R) has been done explicitly up to N = 3 in order to give examples 



for the algorithm defined in (144). How to proceed for N > 3 then should be obvious. 
The reduction rules can be implemented easily in algebraic manipulation programs. 



N=2: 

Necessary for a non-vanishing result is that the external momentum is not light-like. The Gram 
matrix is trivial. One has VV(r) = W u ' 2 = r^r^ 2 /(2r • r). 



I^{r) = — 1 - (nr^V 2 -rV" 12 ) 



(66) 



N=3: 

We use the short-hand notation for the arguments of the tensors and integrals defined in Eq. ( [45| ) 
in the following: R = [ri,r 2 ,0], R[\\ — [r 2 — ri, 0], -R[i] = [r 2 ,Q], etc. The arguments for the tensors 
W*,/Cf , H' 1V are not written explicitly, e.g. it is always understood that W M = W^{R), etc. 



TO 



if 2 (i?) 



w / 3 ™(ii) +^/cr (/ 2 "(ii w ) - / 2 "(i? w )) 

Z=l 

i (>V^ If (fl) + W^ 2 (i?)) 

2 

+ £ (/Cf 1 ?^ 2 + /Cf 2 ?C) + - 

Z=l 

-(/ 1 " > /2-W tt ' ft )Jj +J (fl) 



{A11U2A13} 



Z=l 
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-ir\R[i])-i2 lV2 {R[i])) 



(g/2-nyw 



{^1/^2^3} 



is +2 (R)- (g/2-nyiCt 



{A»lA"2A»3} 



I^(R) (67) 



N=4: 

R = [ri, r2, r3, 0]. The integrals I^il!j L are given in unshifted form. The integrals J^ 1 " 4J ' L are given 
in Eq. @. 



1 



J 4 



T-P1M2A'3A I 4 



W" 1 If (i?) + If (fl)) + ^ \k-{H v 



{M1M2} 



^4-/( i? ) + J 4 lAI2 



(i?) = 



- ( W M1 {If 113 - " 3 } + W 2 {/£ lM3 - J 4 M1M3 } + W" 3 {If M2 - 2 } 
3 



{M1M2M3} 



Z = l 



_ i. ^yyMi |jA*2M3A»4 _ JM2/J3M4 j _|_ yyM2 |jMlM3M4 _ JM1M3M4 j 

}) 



_|_yyA»3 jj^i^/M _ jAHAi2M4| _|_ yyA*4 |jMiA 1 2A'3 _ j^i^2M3 
3 



;=i 



"1 ^"2 ^"3 



{Ml/i2/i3M4} 



JU1V2V3 _|_ JA I 1A'2M3M4 



(68) 



N=5: 

As discussed above, no higher dimensional integrals appear in the reduction of tensor TV-point 
integrals for AT > 4, as long as the external momenta are non-exceptional. They occur only implicitly 
by reduction of 2,3,4-point tensor integrals. Now, R = [rj, r%, T3, T4, 0]. The formulas are valid up 
to 0(e) since terms have been dropped. 



'■(R) 



w^i^R) + J2^ 1 ih- l (R) 

2=1 

, 4 

- (V 1 if (-R) + w^ 2 /f (iz)J + £ 



{^1^2} 



^s-iiR) 



J^ 1 ^ 2 ^ 3 (Jl'\ = _ yy/Mi JM2M3 _|_ yy/J2 jMiW _|_ yy^3 JM1M2 
5 V ■/ 3 ^ ^ 



{M1M2M3} 



TOW 



jA'iM2M3M4 _ _ /'yyA'i jA I 2M3M4 _|_ yyM2 JM1M3M4 _|_ yyM3 JM1M2M4 _|_ yyM4 JA*1A*2M3^ 



+ 



K-'{H Vl H V2 H V3 



{M1A»2M3M4} 



TO"*( fl ) 



18 



1 5 



fJ ' 5 (R) = — (VV^ 1 J^ 12 ^ 3 ^ 4 ^ 5 _|_ _|_ yyM5 jA'iA I 2M3M4^ 



2 

15 



ilT-T Vi {R) 



(69) 



N=6: 

For non-exceptional kinematics any set of four vectors out of {n, ra, ?"3, #"4, ^5} spans 4-dimensional 
Minkowski space. We will express the tensor 6-point functions in terms of scalar 6-point functions 
and pentagon integrals. The latter then can be reduced further by using (|2^) and the reduction 
formulas for N — 5 with the corresponding new argument vectors. 



if 2 (R) 



1=1 



1=1 



It (R) + W C (R)) + £ [W*] bll " 2i ll 6 -i(R) 



1=1 



{M1M2M3} 



TOW 



4 V 



jA'2A 1 3M4 _|_ yyM2 jAHf»sM4 



_l_ yyP3 JM1M2M4 _j_ yyM4 jginzvaj 



{M1M2M3M4} 



TO*(*) 



jMlM2A t 3M4/i5 



jMlA t 2A 1 3A'4MBM6 



5 V 



M2/J3M4AI5 

6 



yy/Js j/'lM2Ai3A I 4 



15 



1=1 



{M1M2M3M4M5} 



I^1T 1 {R) 



yyMl jA I 2/i3M4M5M6 



45 



1=1 



{A I lM2M3A«4M5A I 6} 
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